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Abstract. A bi-Hamiltonian hierarchy of quaternion sohton equations is derived from 
geometric non-stretching flows of curves in the quaternionic projective space HP". The 
derivation adapts the method and results in recent work by one of us on the Hamilton- 
ian structure of non-stretching curve flows in Riemannian symmetric spaces M — G/H hy 
viewing HP" ~ U(n + 1,H)/U(1,H) x U(n,H) ~ Sp(n -1- 1)/Sp(l) x Sp(n) as a symmetric 
space in terms of compact real symplectic groups and quaternion unitary groups. As main 
results, scalar-vector (multi-component) versions of the sine-Gordon (SG) equation and the 
modified Korteweg-de Vries (mKdV) equation are obtained along with their bi-Hamiltonian 
integrability structure consisting of a shared hierarchy of quaternionic symmetries and con- 
servation laws generated by a hereditary recursion operator. The corresponding geometric 
curve flows in HP" are shown to be described by a non-stretching wave map and a mKdV 
analog of a non-stretching Schrodinger map. 



1. Introduction and Summary 

Certain geometric flows of curves in liomogeneous plane and space geometries are well- 
known to encode scalar soliton equations through the induced evolution of geometrical in- 
variants of the curve [HEIEIIIIEIEIEIEIEI [ID]- There has been much recent interest 
in extending such geometrical derivations to mult i- component soliton equations. For ex- 
ample, the two vector versions of the modified Korteweg-de Vries (mKdV) equation known 
from symmetry-integrability classifications [11] have been derived [12] from non-stretching 
curve flows in the homogeneous Riemannian geometries M = SO{N + 1)/S0{N) ~ and 
M = SU{N)/ SO{N), generalizing a similar derivation of the scalar mKdV equation ob- 
tained previously |13j in the case of the standard 2-sphere geometry, ~ S0{3)/ S0{2) ~ 
SU{2) I S0{2). The same approach has also been used to provide a geometric origin for the 
scalar sine-Gordon (SG) equation and its two known vector versions [T3l[T2]. Related work 
on integrable curve flows in Euclidean and spherical Riemannian geometries, and , 
has appeared in [IS [161 [13 [HI [19] . 

The derivation of vector mKdV and SG equations in the geometries M = G/SO{N) is 
based on a moving parallel frame formulation for arclength-parameterized curves 7 which is 
closely analogous to the standard parallel framing of curves in Euclidean geometry [20|, |15l 
[T7] . In this generalization [T2], the Euclidean group is replaced by the respective isometry 
groups G = SO{N + 1) and G = SU{N) whose rotation subgroups SO{N) C G act as 
the gauge group for the frame bundle of M = G/SO{N) (analogously to that of Euclidean 
space). Given a parallel framing, where x is the G-invariant arclength along 7 and T = 
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is the unit tangent vector, the Cartan structure equations for torsion and curvature of the 
Riemannian connections V^, Vt on the 2-dimensional surface of any non-stretching curve 
flow 'y{t,x) in the manifold M = G/SO{N) are seen [12] to geometrically encode a pair 
of compatible Hamiltonian operators. This bi-Hamiltonian structure generates a hierarchy 
of integrable curve flows in which the frame components of the principal normal vector 

= VxT satisfy vector evolution equations related by a hereditary recursion operator. Each 
evolution equation displays invariance under the rotation group SO{N — 1) (in the SO{N) 
gauge group) that acts in the normal space at each point on the curve 7 and preserves the 
form of the connection matrix of the parallel frame, so thus the principal normal components 
have the meaning of geometrical covariants that are determined by 7 up to SO{N — 1) gauge 
freedom (which represents the equivalence group of the framing). The lowest-order flow in 
the respective hierarchies in M = SO{N +1)/ SO{N) ~ and M = SU{N)/SO{N) yields 
two different vector mKdV equations which are SO{N — l)-invariant. In addition, each 
hierarchy also contains a flow that yields a SO{N — l)-invariant vector SG equation arising 
from the kernel of the hereditary recursion operator. The corresponding geometric curve 
flows 7(t, x) in both hierarchies are found to be given by [131 [l2j wave maps and mKdV 
analogs of Schrodinger maps on M = G / SO{N). 

Complex versions of these vector SG equations and vector mKdV equations have been 
derived [21] through applying the same method to geometric non-stretching curve flows in 
the Lie groups G = SO{N + 1),SU{N) themselves viewed as homogeneous Riemannian 
geometries in the standard manner G ~ G x G/diag{G x G). This approach has also yielded 
the two known vector versions of the nonlinear Schrodinger equation (NLS), along with their 
bi-Hamiltonian integrability structure. 

A broad generalization of such results has been obtained in recent work [22] by one of us 
on the Hamiltonian structure of non-stretching flows of curves in homogeneous Riemannian 
geometries (i.e. symmetric spaces) M = G/H, including compact semisimple Lie group 
geometries M = K for G = K x K, H = diag G. For all these geometries, the results 
give a general geometric derivation of group-invariant (multi-component) SG, mKdV, and 
NLS equations along with their bi-Hamiltonian integrability structure consisting of a shared 
hierarchy of symmetries and conservation laws generated by a group-invariant recursion 
operator. 

In the present paper we adapt this derivation to non-stretching curve flows in the quater- 
nion projective geometry M = HP" by means of the symmetric space isomorphisms 

HP" ~ U(n + 1, H)/U(l, H) X U(n, H) ~ Sp(n + 1)/Sp(l) x Sp(n) (1.1) 

given in terms of Hamilton's quaternions H = span(l, i, j, k). Underlying these isomorphisms 
is a basic identiflcation [23] between the quaternion unitary Lie algebra u(n, H) and the 
compact symplectic real Lie algebra sp{n) holding for n > 1. Our main result will be 
to obtain bi-Hamiltonian geometric curve flows in M = HP" yielding an integrable multi- 
component quaternionic version of the mKdV equation and the SG equation. 

A summary of relevant properties of symplectic groups, quaternion unitary groups, and 
the symmetric space structure of HP" is provided in section 2. In section 3, starting from 
the gauge group U(1,H) x U(n, H) of the frame bundle of HP", we extend the moving 
parallel frame formulation to curves in HP" and explain its main geometrical properties. In 
particular, the components of the principal normal vector along any arclength-parameterized 

2 



curve with this framing define geometrical covariants that transform as a quaternion scalar- 
vector pair (or a single quaternion scalar in the case n = 1) with respect to the equivalence 
group U(1,H) X U(n — 1,M.) (which preserves the form of the parallel connection matrix) 
of the framing. In section 4 we next show how the Cartan structure equations for framed 
curve fiows encode two compatible quaternion Hamiltonian operators. We then use these 
operators in section 5 to generate a hierarchy of quaternion Hamiltonian vector fields and 
corresponding quaternion scalar-vector evolution equations with an explicit bi- Hamiltonian 
structure. This yields a quaternion scalar- vector mKdV equation and a quaternion scalar- 
vector SG equation, both of which are unitarily invariant under U(1,H) x U(n — 1,]HI) and 
share a hierarchy of quaternionic symmetries and conservation laws which are generated by a 
U(l, H) X U(n — 1, E[)-invariant recursion operator. In section 6 we derive the corresponding 
geometric SG and mKdV curve fiows in HP" and show that these fiows respectively are 
described by a non-stretching wave map and a non-stretching mKdV map. Finally, we make 
some concluding remarks in section 7. 

2. Algebraic and Geometric Preliminaries 

Here we collect some useful facts about the relevant algebraic and geometric structure of 
HP"" as a quaternion symmetric space. 

The complex symplectic group Sp(n, C) is the group of matrices g in GL(2n, C) that leaves 
invariant the exterior form ziAZn+i + - ■ ■ + ZnAz2n in terms of coordinates (^i, . . . , Z2n) £ C^", 
i.e. 

9 Jn9 — "^n-i Jn — 

Recall that the group of matrices g in GL(2n, C) that leaves invariant the Hermitian form 
z{zi -|- • • • -|- Z2n'z2n IS the complcx unitary group U(2n), i.e. 

9^9 = hn- 

The compact symplectic group is defined by Sp(n) = Sp(n, C) fl U(2n), whose elements have 
the form 

A^A^B'B^In, A^B-B'A^Q. (2.1) 

In the case n = 1, the compact symplectic group Sp(l) has a well-known identification with 
Hamilton's quaternions. The quaternions are defined as an associative non-commutative 
algebra H of dimension 4 over the real numbers, with generators A; that satisfy the 
multiplicative relations 

— — k"^ — —1, ij — —ji — k, jk — —kj — i, ki — —ik — j. 

Quaternions q & M. — span{l,i,j,k) have two main properties: First, every q has a natural 
complex conjugate defined by the conjugation relations 

-i: J = -j: k^-k, T = 1, 

whereby 

Second, each non-zero q has a multiplicative inverse given by q~^ — q/\q\'^, where |g| = {qqY^^ 
is the norm of q. Thus, H is a normed, hypercomplex, division algebra. 
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Proposition 2.1. There is a group isomorphism 

{A B\ f Al+ A2i Bl + B2l\ ~ a ■ . a ■ ■ , ^^ 

^=[-B A) = [-k + k^ l_l■)---9 = A-VB, = A,^A2^^B,,^B2k (2.2) 

under which the compact symplectic group Sp{n) is identified with the quaternionic unitary 
group l]{n,M.) given by all matrices g G GL(n, H) satisfying g^g = In- In particular, Sp(l) 
is identified with the group formed by unit-norm elements in M, 

g = ^ ai - a22j " ^ = + + + ^^'^^ 

where af + + bf + = 1 = \g\'^. 

The group U(n, H) ~ Sp(?T,) acts on H" by right muhiphcation, i.e. for x G and 
g G U(?2, H), X I — > x^ G H'^. There is a natural Hermitian inner product on defined by 

< x,y >:= xy* 

which is invariant under the group action of l]{n,M), i.e. < ^cjiyg >=< x, y > for all 
g G U(n, H). This inner product is the sum of a Euclidean inner product given by 

Re < X, y >= I < X, y > +i < y, X >= i(xy* + yx*) (2.4) 

and an imaginary skew-form given by 

Im < x,y >= I < x,y > -| < y,x >= |(xy* -yx*). (2.5) 

With respect to the Euclidean inner product, there is an orthogonal decomposition of H" 
into real and imaginary parts H" = R" © Q" where 

Q:= {q eM\ q + q = 0} = span{i,j,k) (2.6) 

will denote the set of imaginary quaternions. 

Proposition 2.2. 

1. The Lie algebra g = sp{n) of G = Sp{n) consists of all matrices g G Qi{2n, C) satisfying 

gJn + J„g* = 0, g + g* = 0. 

2. The group isomorphism Sp(n) ~ U(?2, H) induces a Lie algebra isomorphism between the 
symplectic Lie algebra 5p{n) and the quaternionic unitary Lie algebra u(n, H) consisting of 
all matrices g G Qi{n, H) satisfying g + g =0. This isomorphism is explicitly given by 

/ A B\ / Ai + k2i Bl + B2Z\ - A^A -^R -^R^^0 7^ 

S=^-B aJ = ^-Bi + B2Z p^^_p^^,)^^ = ^ + ^3=^^+^2i + B,j+B2k (2.7) 

where B = B*, A + A*^ = 0. In particular, the Lie algebra sp{l) is identified with u(l,E[) = Q. 

3. The Cartan-Killing form on g = sp(n) ~ u{n,M) is given by 

(gi, g2) = tr(ad(f7i)ad((72)) = 2(n + l)tr(gig2) = A{n + l)Re(tr(gig2)) = (gi, ga) (2.8) 

/or gi ^ gi,g2 ^ g2 given by ^2. ?] ). In particular, the negative-definite Cartan-Killing norm 
is explicitly given by 

(g,g)ap(n) = -4(n+l)tr(AA*+BB*) = -4(n+l)tr(AiA* +A2A* +BiB* +B2B*2) = (g,g)u(n,H). 

(2.9) 
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4- For the case n = 1, where g = ( u ^^i^ u ■ + t|2M ^ sp(l) is identified with q = 

a,2i + bij + G u(l, H) = Q, the norm is given by 

- |(g, g)sp(i) = 4 + hl + hl= =< g > . (2.10) 

The Lie group U(n + 1, H) arises geometrically as the isometry group of the quaternionic 
projective space HP". This space consists of the points on the unit sphere in H""*"^ with the 
identification of pairs of points x and xg for every group element g G U(1,H[). Since the 
symmetry group of the Hermitian inner product on H""^^ is G = U(n + 1, H), then the action 
of G modulo the action of U(1,EI) represents the non-trivial isometries of HP". Moreover, 
if we consider the origin in HP" as represented by the point a = (1,0,..., 0) G H"+^, 
then the isotropy subgroup of \J{n + 1,H) leaving invariant this point x = o is given by 
H = U(1,H) X U(n, H). Hence HP" is a homogeneous quaternionic manifold 

-U(l,e)xU(n,H) ^^-^^^ 

on which H* ^ Ad(U(l,H) x U[n,M)) acts as the isotropy group at the origin x = a, i.e. 
H* linearly maps the tangent space TqHP" into itself. It is also clear that the group element 

S := (^~^ G U(r2 + l,e) (2.12) 

yields an involutive automorphism of U(n + 1,H). Consequently, through Proposition 12. H 
HP" has the structure of a symmetric Riemannian space 

-Sp(l)xSp(n)- ^^■^'^^ 

The quaternionic structure of this space takes the form of a triple of linear maps /, J, K on 
T,.HP" having the properties IJ = K,JK = I, KI = J, P = J"^ = K'^ = -id, as canonically 
associated with the action of the isotropy subgroup Ad(Sp(l)) C H* on ToHP". 

At the Lie algebra level, there is a decomposition of = u{n + 1,H) ~ 5p{n + 1) as a 
symmetric Lie algebra constructed in terms of the involutive automorphism cr(g) = SgS. 
The eigenvalues of a are 1 and —1 with corresponding eigenspaces 

f) :=u(l,H)©u(n,H) ~sp(l)©sp(n), a{i)) = I) (2.14) 

and 

m := u(n + l,H)/(u(l,H) ©u(ra,e)) ^ sp(n + l)/(sp(l) © sp(n)), a(m) = -m. (2.15) 

Therefore, with respect to the Cartan-Killing form on g, a induces an orthogonal decompo- 
sition given by the direct sum of vector spaces g = f) © m, with Lie bracket relations 

[f),{)]c{), [f),m]cm, [m,m]cf). (2.16) 

Moreover, the Lie subalgebra P) is identified with the generators of isometries that leave fixed 
the origin in HP" (i.e. yielding the action of H* : T^HP" T^HP"), while the vector space 
m is identified with the generators of isometries that carry the origin a to any point x ^ o 
in HP". 
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Lemma 2.3. 

1. The matrix representation of the vector space m and the Lie subalgebra [) o/g in gl{n+l, H) 
is given by 

(0 m m\ 
-m e m, (2.17) 
-m* / 

/p 0\ /p o\ /o 0\ 
(p,q,h,H):= q h = + q h g = u(l, H) © u(n, H), 
VO -h' H/ VO 0/ \0 -h* H/ 

(2.18) 

m which m, h e H"~-^ are quaternionic vectors, m G EI is a quaternionic scalar, p, q e Q are 
imaginary quaternions, and H G u{n — 1,H) is a quaternionic matrix. The bracket relations 
Ii2.16\) take the form: 

[(m, m), (p, q, h, H)] =(mq — pm— < m, h > , mh — pm + mH) G m (2.19a) 

[(mi, mi), (m2, m2)] =(-mim2 + m2mi - 2Im < mi, m2 > , -mim2 + m2mi , 

— mim2 + m2mi , — m^m2 + rrigmi) G f), (2.19b) 

[(pi, qi, hi. Hi), (p2, q2, ha, H2)] =(piP2 - P2P1 , qiq2 - q2qi - 2Im < hi, h2 > , 

qih2 - qahi + hiH2 - h2Hi , 

h*h2 - hahl + H1H2 - H2H1) G [). (2.19c) 

2. The restriction of Cartan-Killing form on q to m yields a negative- definite inner product 

((mi, mi), (m2, m2)) = -xRe(mim2+ < mi, m2 >) (2.20) 

where 

X = S{n + 2). (2.21) 

3. The rank ofvn. is equal to 1. 

Cartan subspaces of m are defined as a maximal abelian subspace a C m, having the 
property that it is the centrahzer of its elements, o = m fl c(a). It can be shown [21] that 
any two Cartan subspaces are isomorphic to one another under some linear transformation 
in Ad (if) and that the action of Ad (if) on any Cartan subspace a generates m. Since 
m = u(n + l,E[)/(u(l,E[) © u(n, H)) has rank(m) = 1, its Cartan subspaces a are one- 
dimensional. 

A choice of basis element defining a = span(e) is given by 

/O 1 0\ 

e := (1,0) = -1 G m, (2.22) 
\ 0/ 

which can be readily verified to satisfy 

m n c(e) = span(e), m = span(Ad(U(l, H) x U(n, e))e) (2.23) 

where Ad(/i)e = heh^^ for aU h = exp((p, q, h, H)) G U(1,H) x U(n, H) with p, q G Q, 
h G H""^, H G u{n — 1, H). Following the notation used in [22], let the centralizer subspaces 
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of e in m and f) be denoted 

my := c(e) n m, := c(e) n I). (2.24) 

The orthogonal complements (perp spaces) of these subspaces m|| and with respect to the 
Cartan-Killing form will be denoted by rrij. and \)±, where 

m = m|| © m_L, i)^^\\®i)±- 

Their matrix representation is given by 

mil 




(my) := I -my | e m||, (m^,m^):=| m^ | e m^, (2.25) 







m± 


m^ 










'hx 













-hi 



and 

(h||,H||):= I h|| 0|e[)||, (h^,hx):=|0 -h^ I e (2.26) 

\o nj \o -hi / 

where my e M is a real quaternion, h||,mx,hx G Q are imaginary quaternions, mx,hx G 
H""^ arc quaternionic vectors, and Hy G u(n — 1,]HI) is a quaternionic matrix. The corre- 
sponding decomposition of m ~ M"- as a vector space is given by 

my~R, m±~Q©H"-^ (2.27) 

The linear operator ad(e) maps i)± into m_L, and vice versa. In particular, 

ad(e)(h^, h^) = (-2h^, h^) e m^, ad(e)(m^, m^) = (2m^, -m^) e f)^. (2.28) 

Consequently, ad(e)^ is well-defined as a linear mapping of the subspaces l)± and into 
themselves. 

Let H* be the subgroup in the isotropy group H* ~ Ad{H) given by Ad{h)e = e, h G H, 
generated by the Lie subalgebra f)y C f). The subgroup H* can be identified with the group 
U(n — 1, H) X U(l, H) C U(n -I- 1, H) whose matrix representation is given by 

'a 0\ 

a , aeU(l,H), AeU(n-l,H). (2.29) 
.0 Aj 

Here the subgroup U(n — 1, H) acts on m by right multiplication, 

Ad(A)m|| = my, Ad(A)m_L = m_L, Ad{A)m± ^ ni±A, (2.30) 
while U(1,HI) has a non-standard action given by 

Ad(a)m|| = my, Ad{a)m.± — am±a~^, Ad(a)mx = ain.±. (2-31) 
The corresponding group action on f) is given by 

Ad(a)hy = ahya~\ Ad(a)Hy=Hy, Ad{a)h± ^ ah^a"^ , Ad{a)h±^ah±. (2.32) 

Proposition 2.4. The vector subspaces i)± and m± are isomorphic under the linear map 
ad(e), and consequently they each decompose into a direct sum of vector subspaces given by 
irreducible representations of the group H* on which the linear map ad(e)^ is a multiple of 

the identity: 

ad(e)2(mx, m^) = (-4mx, -m^) e m^, eid{ef{h^, h^) = (-4hx, -h^) e f)±. (2.33) 
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The Lie bracket relations on m = m|| © m_L and i) = Q) l)± coming from the structure of 
g as a symmetric Lie algebra fl2.16p consist of 

[m||,m||]C[)||, [m||,f3||] Cmii, [f)||, f^n] C f),,, (2.34) 
[m||,m^] C [)^, [xn\\,i)A_]Cxn±, [t)||, m^] C m^, i)±- (2.35) 

The only Lie brackets with nontrivial decompositions are [mj_,m^], [f)±,t)±], [tn^,f)±]- To 
write out all these brackets explicitly, we introduce the following commutator and anti- 
commutator notations. For a, b G Q and a, b G H""-^, let 

C(a,b) := ab - ba G Q, A(a, b) := ab + ba G M, (2.36a) 

C(a, b) :=< a, b > - < b, a >= ab* - ba* G Q, C(a, b) = a*b - b*a G u(n - 1, H), 

(2.36b) 

A{a, b) :=< a, b > + < b, a >= ab* + ba* G R. (2.36c) 
Lemma 2.5. 

1. The Lie brackets { 2.34 ), f2.35\} are given by 

[(mi||),(m2||)] = 0Gf)||, (2.37a) 

[(mi|),(h||,H||)]=OGm||, (2.37b) 

[(hill, Hill), (h2||,H2||)] = (C(hi||,h2||), [Hi||,H2||]) G [)||, (2.37c) 

[(m||), (m^,m^)] = (2m||m_L , -mym^) G f)_L, (2.37d) 

[(my), (lu, hx)] = (-2m||hx , myh^) G m^, (2.37e) 

[(h||,H||),(mx,m^)] = (C(h||,mx), hymx -mxH||) G m^, (2.37f) 

[(h||,H||),(hx,hx)] = (C(h||,hx), h||hx-hxH||) G [)^. (2.37g) 

2. The remaining Lie brackets are given by the projections 

[(mi_L,mi_L), (m2_L,m2_L)](,|| = (C(mi_L, m2_L) + iC(m2_L, mix), C(m2_L, mix)) G f)||, 

(2.38a) 

[(mix, mix), (m2x, m2j_)]t,^ = (iC(m2x, mix) , mixm2x - m2xmix) G P)x, (2.38b) 
[(hix, hix), (h2x, h2x)]„,| = (C(hix, h2x) - |C(h2x, hix) , C(h2x, hix)) G f)||, (2.38c) 
[(hix, hix), (h2±, h2x)]„^ = (iC(hix, h2x) , h2xhix - hixh2x) G ()x, (2.38d) 
[(mx,mx),(hx,hx)]m|| = ( - A(mx, hx) - ^^(mx, hx)) G my, (2.38e) 
[(mx, mx), (hx, hx)]mx = {lC{h±, m±) , mxhx - hxmx) G mx- (2.38f) 

3. The Cartan-Killing form on i)± ~ mx is given by 

((hi±,hix), (h2x,h2x)) = xRe(hixh2x - hixh^x) = |(A(hix,h2x) - A(hix,h2x)). (2.39) 

An explicit basis for m can be built from the basis {1, i,j, k} for M © Q together with the 
standard orthonormal basis {e/};=i^ „_i of M""^ as follows: 

mil = span(e) ~ M, 

• , . ■ i , -1 (2-40) 

mx = span(ml, m^^, m^) © span(mxz, m^,, m^^^, ml,),=i,..„„_i ~ Q © 



in which 



eA / qer 

tnn = {0,ei)=\ 0, = (0, ge,) = 

-e\ 0/ \e\q 



ml = (g, 0) 




(2.41) 



where q = i,j, k, and / = 1, . . . , n — 1. Note and m^ are, respectively, related to m^; 

and e by the action of U(1,HI) C U(n + defined by Ad(a)e = ml, Ad{a)m±i = ml^ 
where 



1 I , q = i,j,keQ. (2.42) 

,0 In-l^ 

In addition, m±i is itself related to e by the action of the rotation subgroup S0(?2, M) C 
U(n,H) C U(n + 1,H). 

Finally, we record several useful quaternionic identities connected with properties of the 
Cartan-Killing form (12.391) . 

(1) For a G Q, 6, c e H and b, c G H"^^ : 



<b,c> =< c,b>, Re(a <b,c>) = Re(< b,c> a) = -Re(a < c, 6 >), (2.43) 



< b, c > =< c, b >, Re(a < b, c >) = Re(< b, c > a) = -Re(a < c, b >). (2.44) 
(2) For a,b,ceQ: 

Re(a6) = Re(6a), Re(a6c) = Re(6ca) = Re(ca6) = — Re(6ac) = — Re(c6a) = — Re(acfe). 

(2.45) 

3. Frame Formulation of Non-Stretching Curve Flows in HP" 

We consider non-stretching fiows of framed curves in the quaternionic projective space 
M = HP" viewed as a symmetric Riemannian geometry (12.111) . Its Riemannian structure 
will be described in terms of [25j a m- valued linear coframe e and a f)-valued linear connection 
uj whose torsion and curvature 

% := de+[uj,e], ^ := du + l[uj , lu] (3.1) 

are 2-forms with respective values in m and f), given by the following Cartan structure 
equations: 

1 = 0, 9^=-i[e,e]. (3.2) 

Here [■, ■] denotes the Lie bracket on g = u(n -|- 1,H) composed with the wedge product 
on T*M ~ H". This structure together with the (negative-definite) Cartan-Killing form 
determines a Riemannian metric and Riemannian connection (i.e. covariant derivative) on 
the manifold M as follows. For all X, Y in T^M, 

g{X,Y) := -{ex,eY), e\V xY ■= dxCy + [(^x, ey], (3.3) 

where the coframe provides an identification between the tangent space T^M and the vector 
space m = as given by e\X := ex, e\Y := ey G m. The connection is metric compatible. 



V(? = 0, and torsion-free, T = 0, while its curvature is covariantly constant, Vi? = 0. In 
particular, 

e\R{X, Y)Z = [^\ (X A r), ez] = -[[ex, ey], e^], ejT(X, Y) = T(X A F) = 0, (3.4) 

where T(X, F) := V xY - VyX - [X, Y] is the torsion tensor and R{X, Y) := [Vx, Vy] - 
V[x,r] is the curvature tensor. Note the linear coframe and linear connection have gauge 
freedom given by the following transformations 

e — > Kd{h-^)e, uj — > Kd{h'^)uj + h'^ dh (3.5) 

for an arbitrary function h : M ^ H = U(n, H) x U(1,EI) C U(n + 1,EI). These gauge 
transformations comprise a local (a;-dependent) representation of the linear isotropy group 
H* = Ad{H) which defines the gauge group [26] of the frame bundle of M = EIP*^. Both the 
metric tensor g and curvature tensor i? on M = HP" are gauge invariant. 

Let 7(t, x) be a flow of any smooth curve in M = MF"'. We write X = for the tangent 
vector and F = 7^ for the evolution vector at each point x along the curve. Note the flow is 
non-stretching provided that it preserves the U(n+ 1, ]HI)-invariant arclength ds = {'j^ldx, or 
equivalently Vt|7a;| = 0, in which case we can put 

l7x| = 1 (3.6) 

without loss of generality. For smooth flows that are transverse to the curve (such that X 
and Y are linearly independent), 7(t,a;) will describe a smooth two-dimensional surface in 
M. The pullback of the Cartan structure equations fl3.2p to this surface yields a framing of 
the curve evolving under the flow: 

D^et - DtCx + [cu^., ct] - [ut, ej^ = 0, (3.7) 

D^uJt - DtUx + [cux, i^t] = -[ex, et], (3.8) 

with 

ex := e\X = e\jx, et := e\Y = e\jt, (3.9) 

Ux := uj\X = ul'jx, ^^t ■= ^^\y = ^^ht, (3.10) 

where Dx,Dt denote total derivatives with respect to x,t. 

Remarkably, for any non-stretching curve flow, these torsion and curvature equations 
fl3.7p - fl3.10p encode an explicit bi-Hamiltonian structure once a specific choice of frame along 
'j{t,x) is made [22]. A frame consists of a set of orthonormal vectors that span the tangent 
space of M at each point x on the curve. Associated to a frame is the connection matrix 
consisting of the set of frame components of the covariant x-derivative of each frame vector 
along the curve ^7]. Such a framing is obtained from the Lie-algebra components of e* and 
ujx when an orthonormal basis is introduced for m and P) with respect to the Cartan-Killing 
form, where e* is a m-valued linear frame defined to be dual to the coframe e by the condition 
that — (e*, e) = id is the identity map on each tangent space T^M (cf [221I25])- In particular, 
uJx determines e along the curve via the transport equation 

VxC = -ad{ux)e. (3.11) 

Then if {m;};=i . 4„_i is any fixed orthonormal basis for the perp space of span(e2;) in m, 
the set of vectors given by X±i := — (e*,m;) together with X = —{e*,ex) defines a frame 
at each point x along the curve, i.e. span(X_Li, . . . , X_l4„_i, X) = T^M with the subset 
{Xx/}/=i 4„_i being an orthonormal basis for the normal space of the curve. It turns 
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out that the bi-Hamiltonian structure encoded in the resulting frame structure equations 
for 7(t,x) is independent of the choice of a basis, and accordingly the simplest algebraic 
formulation of this encoding arises in terms of the m-valued coframe and ()-valued connection 
variables and (IXTUD . 

We utilize a natural choice of a moving frame defined by the following two properties 
which are a direct algebraic generalization of a parallel moving frame in Euclidean geometry. 

(1) is a constant unit-norm element lying in a fixed Cartan subspace a C m (with 
dim(a) = rank(m) = 1), i.e. Dj-e^^ = DtCx = 0, {e^, = —1, and a = span(ea;) = m|| 
where m|| © m_L = m and (m||, m^) = 0. 

(2) Ux lies in the perp space of the Lie subalgebra [)\\ C f) of the linear isotropy 
group H^^ C H* Ad(U(l,EI) x U(?t,,E[)) that preserves e^, i.e. a.d{i)\\)ex = and 
{ux, f)||) = where f)|| © = f) and f)_L) = 0. 

Such a moving frame for 7(t, x) will be called U(l, H) x U(n, M.)-parallel and its existence 
can be established by constructing a suitable gauge transformation (13.51) on an arbitrary 
moving frame at each point x along the curve Since every Cartan subspace a C m 

is one-dimensional, a U(1,EI) x U(n, H[)-parallel moving frame is unique up to the rigid (a;- 
independent) action of H* on e and Ux- Specifically, given any m-valued linear coframe e and 
f)- valued linear connection matrix uix along 7, we can first find a gauge transformation such 
that h~^exh = e^^ is a constant element in any Cartan subspace a C m, as a consequence 
of the fact m = Ad{H)a. The norm of Cx will satisfy —{ex,ex) = gipfx^lx) = llxl"^ = 1 
because we have chosen an arclength parametrization of the curve. We can then find a gauge 
transformation belonging to the subgroup H* which preserves Cx, so that Dxh+ h~^ujxh = 

ijJx where h{x) E H^^ is given by solving the linear matrix ODE Dxh + roll/i = in terms 

of the decomposition oi tUx = tx^" + relative to Cx- Note the solution will depend on 

an arbitrary initial condition h{xo) G specified at some point x = xq along the curve, 

which represents a rigid gauge freedom (i.e. the equivalence group) in the construction of 
the moving frame. 

Employing the quaternionic matrix notation and algebraic preliminaries from Section 2, 
we choose 

I I ^ ^ ^\ I 

— -1 =: —(1,0) G a 

Vx \ 0) Vx 

(3.12) 



u 



,0 



' 


1 







.0 











— u 


u 


-u* 






(u, u) G [)x ^ Q © H 



n-l 



where u G Q is an imaginary quaternion variable, and u G H"^^ is a quaternionic vector 
variable. (Note the factor 1/^^ in ^x comes from the normalization factor (12.211) of the 
Cartan-Killing inner product on m.) The equivalence group H*^ of the U(1,H) x U[n,M)- 
parallel frame obtained from this choice of and Ux consists of rigid (x-independent) gauge 
transformations (13.51) that preserve the form of the matrices (I3.12p . Specifically, we have 
from flCTl) and (Km . Ad{h-^){1,0) = (1,0) and Ad(/i-i)(u, u) = {u,u) where 

n = ana'^ G Q and u = auA G H"~^ (3.13) 
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for all constant functions (on M) h = (a, A) G H\\ = U(l, H) x U(n — 1, H) with the explicit 
matrix form fl2.29p . 

Hereafter e will denote the m-valued linear coframe along 'y{t,x) as determined (up to 
equivalence) by fl3.12p through the transport condition (I3.1ip . In terms of the dual frame 
e*, the associated framing for T^M will look like 

X := -(e*,e^) = 7^, X^; := -(e*, m^;) ± 7^, 

Xl:=-(e%ml)±7., X^, := -(e*, ml,) ± 7., q = ^,J,keQ ' 

where {m^, m^^, m|^,};=i^ „_i give an orthonormal basis for m_L while {e^} is a unit basis 
for my whence {X, X]^, X^,, X]_Ji 

=i,...,n-i is an orthonormal frame along 7(t,x) such that 
dx = ds is the U(n + 1, H[)-invariant arclength. An explicit choice of such a basis for m is 
shown in (12.411) . The action of Ad(if||) ^ H^^ on this basis will yield all other choices of an 
orthonormal basis for m as given by fl2.30p - (l2.3ip . 

The tangent space of M along '~f{t,x) thereby has an orthogonal decomposition into real 
and imaginary subspaces T^M ^ M" © Q" defined by 

= span(X, X^Ji=i,...,„_i =: ReT^M, = span(Xl, Xl,)^;'^-! =: ImT^M. 

In terms of this structure, the U(1,E[) x U(n, H) -parallel frame (13.140 has the geometrical 
property 

ReVi^X G span(X_|_]^, . . . , X_|_„_j^) and ReV^X^; G span(X), / = l,...,n— 1, (3.15) 

closely resembling a Euclidean parallel frame f2U\ . 

The geometrical meaning of the U(l, H) x U(?t,, ]H[)-parallel connection determined by (13.140 
is seen from looking at the frame components of the principal normal vector 

X:= V,X= (e*,ad(e,)w,) (3.16) 

given by 

^ / 2u -u 

ejX = -ad(e^)u;^. = — I 2u | G m^. (3.17) 



These components u and u are invariantly defined by 7(t, x) up to the rigid (x-independent) 
action of the equivalence group H*^ ~ Ad(U(l, H) x U(n — 1, H)) that preserves the framing at 
each point x. Hence, in geometrical terms, the scalar- vector pair (u, u) describes covariants 
of the curve 7 relative to the group H*, such that each component u G Q, u G H""^ belongs 
to an irreducible representation of this group corresponding to the orthogonal decomposition 
of the vector space ~ Q © HI"^-^. Moreover, x-derivatives of (u, u) describe differential 
covariants of 7 relative to Hp which arise geometrically from the frame components of 
Va;-derivatives of the principal normal vector N. For example, {uxiVLx) corresponds to 

ej Vx-X = -d.d{e^)D^u}^ + ad(t^x)^e^ 

1 / 2ux + 4u^-|u|2 -u^. -3uu\ i?>lK) 

' 2u^ - 4u2 + |up G m^. 




\ u* - 3u*u 

We thus note that the geometric invariants of 7 as defined by Riemannian inner products of 
the tangent vector X = 7a; and its derivatives X = Vxlx-, V^X = V^7x, etc. along the curve 
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7 can be expressed as scalars formed from Cartan-Killing inner products of the covariants 
(u, u) and differential covariants {ux,Ux), {n2x,'^2x), etc.; in particular, 

g{N, N) = -g{X, V^X) = 4|u|2 + \u\^ (3.19) 

g{N,VxN) = -g{X,VlX) = 4Re(u, u,) + Re(u, u,.), (3.20) 

giy.N, V.iV) = g{X, WlX) = 4|u,|2 + |u,f + (4|up + lu^^ + <d\n\M^ + 6Re(uu, u,) 

(3.21) 

comprise all U(?t, + 1, EI)-invariants depending on at most u, u, u^., u^,.. 

Remark 3.1. The set of An — 1 invariants given by {(^(X, V^.X)};=2,...,4n. (and their x- 
derivatives up to differential order An — I) generate the components of the connection matrix 
of a classical Frenet frame [27] determined by jx- 



4. Bi-Hamiltonian Structure 

Let 7(t, x) be any non-stretching curve flow in M = HP" and choose a U(l, H) x U(n, H)- 
parallel framing given by (13.121) . In the flow direction Y = jt, we decompose Ct = h\\ + h± 
relative to Cx, with 

1 / h|| 0\ , 

e mil ~ M, (4.1) 




Gmx^Q©e"-\ (4.2) 

J 

and likewise for cj^ = cc" + vj-^ with 

/w" 

ro" = w" I =: (w",W") G ()|| ~Q©u(n-l,H), (4.3) 



/ w'l 












w'l 





)= 


Vo 





w" 




/w^ 







^ 






-w^ 




u 






J 



U7^=\0 -w^ =:(w^,w^)Gflx^Q©H"~\ (4.4) 

\ -w^* / 

where hy e M is a real scalar variable, h_|_, w", w"*" G Q are imaginary quaternion variables, 
hj_, W"*- G H""-^ are quaternion vector variables, and W" G u(n — 1, H) is a quaternion unitary 
matrix. Here we have inserted a factor l/^/x in which corresponds to the normalization 
factor in Cx, in order to simplify later expressions. 

The frame formulation of the flow 7(t, x) is then given by the general results in [22] as 
follows. 

Lemma 4.1. The Cartan structure equations ( [5*. 7| ), ( f5T5]) for any U(1,H) x J] {n^M) -parallel 
linear coframe e and linear connection uj pulled back to the two-dimensional surface 7(t,a:) 
take the form of a flow on uOx ='■ u{t, x) given by 

Ut = D^w^ + [u, ctj"] + [u, w^]± + ad{ex)h±, (4.5) 
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where 

zu^ = ~ad{e^)-\D^h^ + [u, h^] + [u, h±]±), (4.6) 
= roll = -D;\u,w% (4.7) 

are given in terms of h±, with D^^ denoting the formal inverse of the total x -derivative 
operator . 

Note that h± geometrically corresponds to the normal part of the flow vector Y , 

where Y± is the orthogonal projection of Y relative to the tangent vector X along the curve. 
Similarly, the tangential part of the flow vector Y corresponds to /ly, which is related to 
h±_ through (14. 7p as a geometrical consequence of the non-stretching property of the flow, 
Vt|X| = 0, and the torsion-free property of the framing, VjX = ■ 

We thus emphasize that the formulation stated in Lemma WA\ applies to all non- stretching 
curve flows 7(t,a;) in HP", with the flow being determined by specifying Y_[_ freely as a 
function of t at each point x along the curve. 

We now proceed with presenting the bi-Hamiltonian operators encoded in the Cartan 
structure equations (H^ . fHT]) . Write 

/2h^ \ 
■= ad(e,)/ix = — -2h^ -h^ := (h^,h^) G f)x (4.8) 
V h* / 

which belongs to the same space as m = (u, u), where 

Vxh^ = 2h^ G Q, = -h^ G H"-^ (4.9) 

Then the flow equation (14. 5p can be written in operator form [22j 

ut = n{w^) + h^, w^ = J{h^), (4.10) 

where 

n = JC\^^, J = -ad(eJ-^/CU^ad(eJ-^ (4.11) 
are linear operators which act on valued functions and are invariant under H*^, as deflned 
in terms of the operator 

}C:=D^ + [u,-U-[u,D-'[u,-]\\]. (4.12) 

To display these operators explicitly, we flrst write out the flow on the scalar and vector 
variables u = (u, u), yielding 

= D^w^ + C(u, w") + iC(u, w^) + h^, (4.13) 

Ut = D^w^ - wllu + uWll + w-^u - uw^ + h-^, (4.14) 

where 

w" = -D~^ {C{m, w^) - iC(u, w^)) , W" = D-^C{vL, w^), (4.15) 

w^ =x(^I^.-h^ + ^C^(u,h^)) +h||U, w^ = x(l)xh^ + ih^u + uh^) +hnu, (4.16) 

hy = -xD-' (iA(u, h^) - lA(u, h^)) . (4.17) 
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Next we introduce some operator notations. For v G Q, v G H" V G u{n — 1, H), let 

av := C(u, v) G Q, := iC(u, v) G Q, ^ 

CuV := C(u,v) G u(n- l,e), := iA(u,v) G M 

denoting commutator and anti-commutator operators, and let 

RuV := vu G H"-\ LuV := uv G W'-\ L^V := uV G W'-^ (4.19) 

denoting right and left multiplication. Last it is convenient to scale the variables 

h^ ^ -h^, ^ -h^ (4.20) 

X X 

in order to absorb all x factors in ( 14.16p -( l4.17p . 



Theorem 4.2. The scaled flow equations Iji4-13^ , fi-14 ), ^i4^'2^ ) for the quaternion variables 



u{t,x) G Q and u{t,x) G H" ^ have the operator form 



where 



and 



Ti: = I I (4.22) 



= I I (4.23) 

are compatible Hamiltonian cosymplectic and symplectic operators on the x-jet space o/(u, u). 

The proof of Theorem 14.21 follows directly from general results proven in [22] on the 
Hamiltonian structure of non-stretching curve flows in symmetric spaces, as applied to HP" ~ 
Sp(n + 1)/Sp(l) X Sp(ri). The work in [22] also develops the basic theory and properties 
of bi-Hamiltonian operators for Lie-algebra valued flow variables, generalizing the standard 
treatment for scalar variables given in [2B] (see also [25]). 

To explain the definition and properties of Hamiltonian operators in the setting of quater- 
nions, we start by defining variational derivatives on the x-jet space J°° := (x,u, u, u, 
nx,Ux,u.rc, ■ ■ ■) of the quaternionic flow variables u{t,x) and u{t,x). For any real-valued 
functional = J H{x, u, u, u, u^;, u^,., u^., . . .)dx, its variational derivatives with respect to u 
and u are defined in terms of Frechet derivatives of H{x, u, u, u, u^,, u^,, u^,, . . .) by 

S,.^ ^ /p.(h- . 3/9um. ^: /Re(h^(.«/^u)).. ^ Re/ < h^ -SS,/Su > ... (4.24) 

5^^f, = J priReih^ ■d/du))Hdx=: j Re(h^((5i3/5u)') = Re j <\i-^ ,5S)/5vl> dx 

(4.25) 

modulo total x-derivatives, holding for all imaginary quaternion functions h-*- and all quater- 
nionic vector functions h"*", where h-'--(9/(9u and Re (h"'"-9/9u) = |(h"'" -S/Su + h -d/du) are 
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corresponding vector fields prolonged to J°°. (Here the dot denotes summation over quater- 
nion components.) We note these definitions of 6S)/6n and 6Sj/6u involve a reordering of 
products of quaternion variables (as given by the middle equalities of fl4.24p and fl4.25p ). for 
which we use the quaternionic multiplication and conjugation identities fl2.43l) - fl2.45l) . 

Now the property stated in Theorem 14.21 that the operator Ti. is cosymplectic means it 
defines an associated Poisson bracket 

«.} :^ «e / < (-/^«;/^") . H {-f;/,^) > (4.26) 

which is skew-symmetric and obeys the Jacobi identity, for all real- valued functionals S)i,S)2- 
A counterpart of the Poisson bracket is the symplectic 2-form defined in terms of the operator 
J by 

u,(x,x,u:=Re/<(^;;;).j(^j;;)>d. (4,27) 

where Xi,X2 are vector fields X = h.-^ ■ d/dn + Re(h"'" ■ d/du) associated to pairs of an 
imaginary quaternion function h-*- and a quaternion vector function h"*". The property stated 
in Theorem 14.21 that is symplectic corresponds to a; being skew-symmetric and closed. In 
particular, closure means that 

= pr(Xi)a;(X2,X3) + cyclic 

= Re y < (^^^ , pr(h^ ■ d/du + Re{hi ■ d/du)) J >dx + cycfic ^^'^^^ 

holds modulo total x-derivatives for all vector fields Xi, X2, X3. Compatibility of these oper- 
ators H and is the statement that every linear combination ciTC + C2Jl~^ is a cosymplectic 
Hamiltonian operator, or equivalently that CiH~^ + C2J is a symplectic operator, where 7i~^ 
and denote formal inverse operators defined on the x-jet space J°°. 

5. Bi-Hamiltonian Hierarchies of Soliton Equations 

Composition of the compatible Hamiltonian operators (14.221) and fl4.23p yields a recursion 
operator 

(^:; 

given by 

^12 = + C^Dx + CuLu + ^DxA^D^ ^A^^ + C^D^ ^C^Dx — jC^D^ ^C*uC*u 

^21 = \DxRvL + \RuDx — \LnRvL + \RuDx ^CuDx — ^DxRuDx^A^i — i_RuAu-D^Mu 

— ^RuD,j. ^C^R^ + ^L^D^ ^Cu-Ru + ^L^^R^D^ ^^uj 

7^22 = D1 + DxL^ — L^Dx — L^iL^ + ^RuCu + DxRuDx^A^ — R^D^^CuDx + L^D^^C^Dx 

+ \R\iDx C'uCu + \R\iA^xD ,j. — R^Dj. CuL^ + L^Dj. Cu-^u ~ L^R^D^ A^- 

Each of these operators 7i, J7, TZ displays obvious symmetry invariance under translations 
in X. As a consequence, from general results due to Magri [301 EI], the recursion operator 
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will generate a hierarchy of commuting Hamiltonian vector fields with respect to the Poisson 
bracket, starting from the evolutionary form of the x-translation vector field d/dx. Moreover, 
the adjoint recursion operator TZ* will generate an involutive hierarchy of variational covector 
fields, arising from the canonical pairing provided by the symplectic 2-form. This leads to 
the following results 



Theorem 5.1. The pairs of quaternionic scalar-vector functions h^^ G Q, hj-^^ G H" given 
by 

h|)^-^'(u:)' 0.1.2.-. (5-2) 

yield a commuting hierarchy of Hamiltonian vector fields h^^ ■ d/dn + Re(h([) ■ d/du). In 
particular, there exists corresponding Hamiltonian functional S^^''^ such that 

for all functionals (J5 on J°° . Explicit expressions for the Hamiltonians if'-'-' G M, / = 
0, 1, 2, . . are given by 



(5.3) 



H^'^ ■■ = Y^i^^'^<< ^'^) > + < "'hi) >) 

= -^D'\-A(nX)) + M^M))) = 

2 + 4/ ^ ^ ' ^ ' 1 + 21 II 

whose variational derivatives 

wj) := 5H^^^/5u = -5H^^^/5n G Q, w^^ := SH^^^/Su G H""^ 
yield an associated hierarchy of involutive covector fields w^^ ■ du + Re(w^^ ■ du) given by 



w; 



(0 ^ —T>*ir^ 



These variational covector fields are dual to the Hamiltonian vector fields via the pairing 

(wj) ■ dn + Re(wJ) ■ du))\X = a;(X, hfi^ ■ d/dn + Re(hJ) • d/du)) (5.5) 
holding for all vector fields X in evolutionary form on J°°. 

Note both hierarchies (15. 2p and (15.41) possess the mKdV scaling symmetry x Xx, 

(u, u) — i> (A~^u, A~^u), with (15.21) and (15.31) having the scaling weight 2 + 2/ and (15.41) 
having the scaling weight 1 + 2/. 

We can now state our main result. 

Theorem 5.2. The flow equations ^4-21 ) on the imaginary scalar quaternion variable u(t,a;) 
and the n—l-component quaternion vector variable u(t, x) yield a hierarchy of bi- Hamiltonian 
evolution equations 
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/ = 1,2,..., (5.6) 



called the {+l)-flow, with the Hamiltonian structure 

Hi) 



where S := TZTi. Each of these multi- component quaternionic evolution equations Ii5.6\) is 
invariant under the group H^^ ~ Ad(U(l,H) x U(n — 1,IHI)) given by rigid (x-independent) 
transformations 1^3. 13\) on the pair of scalar-vector variables (u, u). 

5.1. mKdV flow. The +1 flow in the hierarchy (15. 6p comes from 



no) 



producing a coupled system of quaternionic scalar-vector mKdV equations 



(5.9) 



- X = iUsx - fu^u^ + |C(u, C(u, u^)) + |C(u, ua^) = hf^^, (5.10) 

ut - = + |(|up - u^ + u^)u^ + f (2u|up - A{n, u^) + - C(u, u^))u = h^-^), 

(5.11) 

where C{-, ■) and A{-, ■) respectively denote the commutators and anticommutators defined 
in fl2:36al) - (l236cl) . 

This is an integrable system in the following sense. We first note that the convective 
terms u^; and u^, in both equations (I5.10l) - (l5.1ip can be removed by a Galilean transformation 
t ^ t,x x + x^^t- The resulting evolution equations then have the explicit bi-Hamiltonian 
structure given by the / = 1 case of fl5.7l) - fl5.8p : 



uA _ ^ f6HW/6n\ _ . f5H(^^/6n 
Ut) ~ \6H^'^/6nJ ~ ^ \6H('^/6u 



lu3x - fu^U:, + f C(u, C{U, Ux)) + |C(U, U2:,.) 

Usx + |(|up - u^ + nx)ux + |(2u|up - A{u, u^) - C{u, u^.) + U2a-)u 



(5.12) 



where 

= -lu2+ |u|2 and i7W = |u^.-l|u,|2-|A(u,C(u,u,)) + |(u2-|un2 (5.13) 
are the Hamiltonians. 

In addition to its explicit symmetry with respect to space translations x ^ x + e and 
mKdV scalings x — > Xx, t — X^t, u — A~^u, u — >■ A~^u, the quaternion scalar-vector mKdV 
system f l5.12p possesses higher symmetries given by each higher-order flow in the hierarchy, 
i.e. X = h^^ -8/811 + Re(h(^) ■ 8/8u) generates an infinitesimal symmetry of the coupled 
system f l5.12p for all I = 2,3,.... Moreover, each Hamiltonian in the hierarchy yields a 
conservation law for this system (15.120 . i.e. H^^^dx = for / = 0, 1, 2, . . . holds for all 

solutions n{t,x),u{t,x) that have sufficiently fast decay as x ^ ±00. 
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5.2. SG flow. Apart from the +1, +2, . . . flows in the mKdV hierarchy (15. 6p . the recursion 
operator TZ = TiJ also yields a flow deflned by 



^=[l^]=J{l^]- (5-14) 



This will be called the —1 flow f22]. The resulting flow equations (14.211) have the form 

"'^-X-'C;') (5.15) 



VLt J \\l 

with 

D^h^ = -C{u, h^) - 4h,|U, D^h^ = -|h^u - uh^ - hyu (5.16) 

and 

D^h|| = -iA(u, h^) - iA(u, h^). (5.17) 

These equations (I5.16p - (l5.17p possess the conservation law 

D,(hJ + i|h^|2 + |h^n = 0. (5.18) 

Hence, after a conformal scaling of t, which induces a corresponding scaling of the variables 
h||, h"*-, h"*" by a function of t, we get 

hj + i|h^|2 + |h^|2 = c = const. (5.19) 

yielding the relation 

h|| = ±^c- i|h^|2 - |h^|2. (5.20) 

Substitution of (15.201) and (15.151) into (I5.16P then gives a coupled hyperbolic system of quater- 
nionic scalar-vector SG equations 



ut. = t4a/(1 - iWt\^ - |ui|2) u - C(u, Ui), (5.21) 



utx = T\l (1 - i|utp - |ut|2) u - iufU - uui, (5.22) 

in which we have put c = without loss of generality and scaled out a factor x by means 
of the transformation x X^^^- Here C{-, ■) denotes the commutator defined in fl2.36bl) . 

This system (I5.2ip - (l5.22p has explicit symmetry under separate time translations t t + e 
and space translations a; — > x + e as well as scalings t — > X~^t, u A~^u, u —>■ A^^u. It also 
possesses higher symmetries X = h^-j -8/8x1 + Re(h(^) ■ 8/8u) for / = 2, 3, . . ., given by each 
higher-order flow (15.20 in the quaternionic mKdV hierarchy (15. 2p . along with conservation 
laws I H^^Ux = for / = 0, 1, 2, . . ., given by the corresponding Hamiltonians (15.31) in 
the same hierarchy, as verifled through (I5.15p - (l5.17p . 

5.3. Reductions and Soliton equations in HP^. We now make some remarks on reduc- 
tions of the quaternionic mKdV system (I5.12p and quaternionic SG system (I5.2ip . 

Firstly, if we consider a reduction to a vector system by putting u = 0, then the quater- 
nion vector inner product terms C(u, U2a;) and C(u, Ut) in the respective equations (15.120 
and (I5.2ip for the quaternion scalar variable u need to vanish identically. Such algebraic 
constraints cannot be satisfied unless we restrict the vector variable u to be commutative 
by taking its components to belong to a real or complex subalgebra of H. Therefore, both 
the quaternion mKdV system (I5.12p and quaternion SG system (I5.2ip have no consistent 
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non-commutative vector reduction. The underlying reason why the quaternion vector vari- 
able u is unavoidably coupled to the quaternion scalar variable u can be understood from 
the nontrivial Lie bracket relation [(0,hij_), (0, h2±)]|| = (|C(hij_, h2_L), O) for the space l)± 
in which u hes (cf. ( HMd^ and flSTT^ l 



In contrast, we can get a consistent scalar reduction by putting u = in the systems 
f l5.12p and fl5.2ip . This reduction has a natural geometric meaning if non-stretching curve 
flows are considered in a submanifold HP^ C HP"' or in HP^ itself when n = 1. 

Thereby we obtain a scalar non- commutative mKdV equation 

Ui = |u3^ - fu^Uj. (5.23) 
and a scalar non-commutative SG equation 

nt, = 4u^l + iu2 (5.24) 

in which n{t, x) is an imaginary quaternion variable. We note the underlying Hamiltonian 
structure in this case is given by Theorems 15.11 and 15.21 with u = and h"*" = 0, where the 
Hamiltonian operators take form 

n = D,- C^D-'C^, J=\a,- \A^D-'A^ (5.25) 

(which are the first diagonal entries of f l4.22p and f l4.23p . respectively). 

6. Geometric Curve Flows 

The bi-Hamiltonian flows given in Theorems 15.11 and 15 . 21 have a direct geometrical formula- 
tion in terms of the U(l, H) x U(n, H[)-parallel frame variables (13. 9p . (I3.10p . (I3.12p describing 
a non-stretching flow of a curve 7 in M = HP". To begin, we recall that the quaternionic 
components of the connection matrix uj^. = (u, u) = u describe covariants of 7 relative to the 
equivalence group of the framing. This means that u is invariantly determined by 7 up to 
the action of the group 

i/,]- ^ Ad(U(l,H) X U(n- 1,H)) (6.1) 

as given by the rigid (x- independent) transformations fl3.13p . The evolution of u under the 
bi-Hamiltonian flows fl5.6p - fl5.8p is given by 

Ut = hfi^,^+X-'hf,^ (6.2) 

with 

hfi^ = n'u^ (6.3) 

for alH = 0, 1, 2, . . ., where TZ is the recursion operator defined in terms of u through (14.110 
and (14.120 . Each evolution equation (16.20 is invariant with respect to the transformations 
(I3.13P and belongs to the general class of flows in which 

h-^{x, u, Ux, Uxx, ■■■) = M{a~^)h^{x, M{a)u, Ad(a)M^, Ad(a)u^a;, . . .), 

a G U(l,e) X U(n- i,e) ^ ' ' 

is an equivariant function of the invariant arclength a; of 7 and the (differential) covariants 
u, Ux, Uxx, ... of 7 relative to the group (16. ip . 

Any flow specified by the U(1,EI) x U[n — 1, EI)-equivariant class (16.40 yields a non- 
stretching curve flow in M = HP" via the geometric relations 

= x'^d.d{ex)-^h^ = ejFx, h = -D'^iu, = e\Y\\ (6.5) 
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taking into account (14. 8p and fl4.20p . where Y± and Y\\ are the normal and tangential projec- 
tions of y = 74 relative to the tangent vector X = 7^. along 7. In particular, the evolution 
vector of the curve is given by 

Y = -{e*,h^ + h\\) = x-'{e*, y{h^)) (6.6) 

in terms of the operator 

y := D-'[u, ad(e,.)-'-]|| - ad(e,.)-' (6.7) 

where e* is the linear frame dual to the linear coframe e along 7, with = e\X. In this 
relation (16. 6p . Cx is preserved under the action of the equivalence group (16. ip . while up to 
equivalence, both e* and e are determined by uj^ through the transport equation (13.110 along 
7. Hence we obtain the following geometric characterization of the class of flows (16.40 . 

Proposition 6.1. Every U(1,H) x U[n — l,M.)-equivariant flow ( [g.^[ ) determines a non- 
stretching curve flow 'yit,x) in HP" satisfying a \J{n + 1,IS) -invariant evolution equation 

7t = y(x,7x,V.7x,V^7x,...) (6.8) 
through Ii6. 6\) -i 6.1). where x is the U(n + 1,M.) -invariant arclength along the curve. 



Invariance of an evolution equation (16. Sp can be shown to imply that the function Y is 
constructed using only the metric g given by (13. 3p and the tensor ad^ defined as follows: 

e\adl{X)Z := ad{exfez (6.9) 

for all X, Z in T^M. Note both g and ad^ are gauge-invariant under changes of frame (13.50 
and thus are well-defined as geometrical structures on the manifold M = HP". 

We can now formulate the bi-Hamiltonian flows from theorem 15.11 and 15.21 in strictly 
geometrical terms. 

Theorem 6.2. The hierarchy of bi-Hamiltonian quaternionic flows / l5.6]) -/ r57^) correspond 
to non-stretching geometric curve flows in M = HP" given by evolution equations of the 
form 

7i = ^(0(7^^, V^7x, V^7x, . . .), |7x| = 1, (/ = 1,2,...) 
where each equation is invariant with respect to the isometry group U(n + 1,H) of M = HP" 
and preserves the invariant arclength x. 

To write out these geometric curve flow equations explicitly, it is useful to introduce the 
linear map 

:= -ad^(7,) (6.10) 

which is determined by the curve 7. In a U(1,H) x \][n — 1, H)-parallel frame, this map 
corresponds to —ad^CxY under which the vector space m = u(n + l,H)/(u(l,H) © u(n,H)) 
decomposes into a direct sum of eigenspaces mj| = M, m]^ = Q, = H"^^ with respective 
eigenvalues 0, 4/x, 1/x, where m|| is the centralizer space of in m, and m_L = © is 
the perp space of m|| . Since the linear coframe e provides an identiflcation between m and 
TxM, there is a corresponding decomposition of the tangent spaces T^M along 7 given by 

T^M = (T^M)i, © (T,M)^, (T^M) II = span(7,), 

(T^M)x = span(7,)^ = (T,M)1 © (T^M)l 
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where {T^M)\\ ~ M, (T^M)^ ^ Q, {T^M)l ~ H"-^ are the eigenspaces of with eigenvalues 
0, 4/x, 

Hereafter, we will write the tangent vector of 7 as T = G (TyM)|| and the principal 
normal vector along 7 as = Vx'Jx ^ (T-yM)^; iV** and will denote the projections of N 
into (T^M)^ and (T^M)^ respectively. 

6.1. mKdV Curve Flow. From the flow (15.91) . after undoing the scaling (14.201) . we have 
y/xi^± = |u^, yxh_L = -u^. through I^Mi, and hy = -^u^ + ||u|2 through (KTl} . The 
frame variables (14. ip and (14. 2 p thereby yield 

. / iu, 0\ 1/00 
(eOl = ^ K h (eOl = ^ , (6.11) 
VX\0 00/ VX\utOO/ 

(et)ii = (-iu' + i|up)e,. (6.12) 

These quaternionic matrices can be expressed in terms of the geometrical vectors T, A^, 
VxN =: N' as follows. We first note from fIXTTD and fICT]) : 

, /O 2u 0\ , /O -u\ 

(e;v)l = — 2u =eJiV^ (e^v)! = — = eJiV^ (e;v)|| = 0, 
\ 0/ \u* / 

(6.13) 



and 



^ / 2u^ 0' 

(e^Oll = (4u' - \u\'')ex, (e,v')l = — 2u, I = eJ(V.iV)l, (6.14) 



0, 



— — 3uu 



{eN')l = —\ I =eJ(V,.iV)l. (6.15) 

\u* -3u*u 



In addition, we note: 



-6uu^ 



(ad(eAr)'e,)± = — ^ | = ej (ad^(iV)T)l. (6.16) 

\-6u*u 

Since |u|2 = ^(A^^, A^^) = jA^^p and -4u2 = g{N',N') = jA^^p, then from (16121) and flXTOD 
we find 

(eOll = + ||Ar''P)e. = ej ((||Arf + ||Ar^p)T). (6.17) 
Next, comparing fleUD to fleTTij) . fl6J[5D . fl616|) . we obtain 

(et)l = i(e^01 = eJ(iV.Ar)l (6.18) 

and 

(et)l = (e^vOl - ix(ad(e^)2e.)l = eJ ((V.AT)! - ix(ad^.(Ar)T)l) . (6.19) 

Thus, from (16.170 . (16.181) . (16.190 combined with et = ej7t, we derive the following geomet- 
rical evolution equation 

7i = IN'I + N'l - ix(ad^.(Ar)r)l + (||Arf + '-\N^\')T (6.20) 
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in terms of T = N = V^T, N' = VxN. We can write this evolution equation in an 
equivalent form without the s and v projections, by relating (16.111) to a.d{ex)~^ applied to 
f l6.13p . (16.141) . (I6.15p . This leads to the geometrical evolution equation 

X7t = V,(A'-V,7x) - (ad2(A'-V,7x)7x)^ - 3(ad^(A'- V,7x)7x) = 1, (6.21) 
or equivalently 

7i = X-'{x-'^h. - |ad^(V.7x)7x)^ - h^'ai^^^'^xl., V.7x)7x, b.\ = 1, (6.22) 

called the non-stretching mKdV map on M = HP"". Each equation (16.201) . (I6.2ip . (I6.22p is 
invariant under the isometry group U(n + 1, H) of M = HP". 

6.2. SG Curve Flow. The flow given by (I5.14p has w^ = w-*- = which implies w" = 
W" = from (14. 5p . Hence the frame variables (14. 3 p and (14.41) yield 

(c^i)± = (^t)ll =0, (6.23) 

so thus the connection matrix in the flow direction vanishes. This can be expressed geomet- 
rically in terms of the tangent vector T = 7^. as follows. By applying a.(l{ex) to Ut = 0, we 
get 

= ad(e,)cut = -[ut, ex] = ej (V^T) (6.24) 
through using DtCx = 0. As a result, we obtain the geometrical evolution equation 

= VtT (6.25) 

or equivalently 

= Vt7x = V.7i, |7x| = 1, (6.26) 
which is called the non-stretching wave map on M = HP". In addition to satisfying the 
non-stretching property Vt|7x| = 0, this equation (I6.26P possesses the conservation law 
Va;|7j| = 0, corresponding to (15.180 . Thus, up to a conformal scaling of t, the evolution 
given by (16.261) describes a flow with unit speed, |7t| = 1. 

The wave map equation (I6.26P and its conservation laws are invariant under the isometry 
group U(n + 1, H) of M = HP". 

6.3. Reductions and Curve Flows in HP^. We remark that the mKdV and SG curve 
flows given by the geometric map equations (I6.22p and (16.260 have a consistent reduction 
such that 7(t,x) is a map into a submanifold HP"*^ C HP" or into HP^ itself when n = 1. In 
the resulting curve flows, the components of the principal normal vector = Vx'~fx along 
the curve 7 in a U(1,H) x U(n — 1, H)-parallel framing satisfy the scalar non-commutative 
mKdV and SG equations (15.231) and (15.241) respectively. 

7. Concluding Remarks 

Our derivation of quaternionic soliton equations (15.60 and their Hamiltonian structure 
(I5.7p - (l5.8p from geometric curve flows (16.80 in HP" can be reformulated entirely as an al- 
gebraic method at the level of the symmetric Lie algebra structure (I2.14p - (l2.15p associated 
with the compact real symplectic group Sp(n) and the quaternionic unitary group U(n, H). 
Specifically, as shown in [22j , the Cartan structure equations (I3.70 - (I3.100 for a framed curve 
flow in any symmetric space M = G/H arise directly from the zero- curvature equation sat- 
isfied by the left- invariant g-valued Maurer-Cartan 1-form ujg on the Lie group G viewed 
as a principal if-bundle over the manifold M. In this setting the puUback of uq by a local 
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section ij) : M G yields the m-valued linear coframe e and ()-valued linear connection 
1-form u; on M via e + uo = iIj*{ujg), where a change in ip ^ ip = iph corresponds to a local 
gauge transformation fl3.5p on e and u. 

Such a zero- curvature approach has been used in recent work [32] to derive a quaternionic 
mKdV equation for the variable u = = i^JTx, based on the choice of a quaternionic 
connection matrix given by 

/O 0\ 
i^x = I — u u I 
\0 -u* 0/ 

in terms of the imaginary scalar quaternion u and the n — 1-component vector quaternion u. 
This differs compared to our choice given by a parallel framing fl3.12p and leads to a more 
complicated bi-Hamiltonian structure p2] 

where 7i and J are a compatible pair of Hamiltonian cosymplectic and symplectic operators 
and A/" is a non-trivial Nijenhuis operator [291 El], producing a different form than flS.lOp - 
(15. lip for the quaternionic scalar-vector mKdV equation arising from the flow defined by 
/i"*- = Uj.. However, it was not shown in [32] whether the above choice for the connection 
matrix can be achieved by a gauge transformation starting from an arbitrary form oi u = 

Building on work by one of us [33], we will show in a forthcoming paper [31] that these 
bi-Hamiltonian operators are in fact related to our operators f l4.22p and f l4.23p by a Backlund 
transformation that can be interpreted as a Hasimoto gauge transformation corresponding to 
a change of framing. We will also derive Lax pairs for these quaternionic mKdV equations, as 
well as the quaternionic SG equations (I5.2ip - (l5.22p . by means of the zero-curvature equations 

dSlD-dsloD. 

In addition, by expanding the imaginary scalar quaternion u and the (n — 1 component) 
vector quaternion u in a quaternionic basis {1, i, j, fc}, we plan to compare the resulting cou- 
pled integrable systems of 4n— 1 ordinary (real- valued) scalar variables to other constructions 
of integrable multi-component mKdV and KdV type systems known in the literature. 

Finally, in another direction, we also plan to explore the possibility of deriving octonion 
soliton equations from bi-Hamiltonian geometric curve flows in the octonion projective plane 
OP^ by adapting the zero- curvature/ moving frame method used in the present paper. 
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